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ABSTRACT 

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are < —21 mag can be 
explained by the ®®Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta- 
circumstellar medium (CSM) interaction model. Recently, iPTFlSehe challenges these energy-source 
models, because the spectral analysis shows that ~ 2.5Mq of ®®Ni have been synthesized but are inad¬ 
equate to power the peak bolometric emission of iPTFlSehe, while the rebrightening of the late-time 
light-curve (LC) and the Hct emission lines indicate that the ejecta-CSM interaction must play a key 
role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar 
together with some amount (< 2.5Mq) of ^®Ni may power the early LC of iPTFlSehe while the late¬ 
time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we 
suggest that iPTFlSehe is a genuine core-collapse supernova rather than a pulsational pair-instability 
supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on 
their explosion and energy-source mechanisms. 

Subject headings: stars: magnetars - supernovae: general ~ supernovae: individual (iPTFlSehe) 


1. INTRODUCTION 

Supernovae (SNe) are exceptionally brilliant stellar ex¬ 
plosions in the Universe. Over the last decade, dozens of 
superluminous SNe (SLSNe) whose peak magnitudes are 
< —21 mag (Gal-Yam 2012) have attracted intense inter¬ 
est from both observational and theoretical astronomers. 
SLSNe are the rarest class of supernovae discovered so far 
and have been classified into types 1 (hydrogen-poor) and 
II (hydrogen-rich). Comparing the results of Cooke et al. 
(2012), Quimby et al. (201S), and Taylor et al. (2014), 
the ratio of the explosion rates of SLSNe to core-collapse 
SNe (CCSNe) is ~ (1 - 4) x lO'^. 

Discriminating among the energy-sources powering 
the SLSNe is rather tricky. However, these energy 
sources should leave their imprints on the light curves 
(LCs) and the spectra of these SLSNe. Only very 
few SLSNe (e.g., SN 2007bi) can be regarded as 
the “pair instability SNe (PlSNe)” (Barkat et al. 
1967; Rakavy & Shaviv 1967; Heger & Woosley 2002; 
Heger et al. 200S) and explained by the ®®Ni-powered 
model (Gal-Yam et al. 2009)^, while the majority of 
SLSNe cannot be solely powered by ^®Ni (Quimby et al. 
2011; Inserra et al. 201S) but can be explained by 
the ejecta-circumstellar medium (CSM) interaction 
model (Chevalier & Irwin 2011; Ginzburg & Balberg 
2012) or magnetar-powered (ultra-highly magnetized 
pulsar) model (Kasen & Bildsten 2010; Woosley 2010; 
Dessart et al. 2012; Inserra et al. 201S; Wang et al. 

^ However, Dessart et al. (2012) and Nicholl et al. (2013) disfa¬ 
vored the ®®Ni-powered model in explaining SN 2007bi and other 
similar events. 


2015a). Because of the absence of emission lines in¬ 
dicative of the ejecta-CSM interaction, the magnetar 
model is preferred in explaining almost all SLSNe-1 (e.g. 
Inserra et al. 201S; Nicholl et al. 201S; Howell et al. 
201S; McCrum et al. 2014; Vreeswijk et al. 2014; 
Nicholl et al. 2014; Wang et al. 2015a; Metzger et al. 
2015; Dai et al. 2016; Kasen et al. 2016; Wang et al. 
2016). Meanwhile, the ejecta-CSM interaction model 
works well in explaining SLSNe-H, especially superlu¬ 
minous SNe Hn (Smith & McCray 2007; Moriya et al. 
201Sa). 

In principle, ®®Ni should be taken into account in light- 
curve modeling for all SLSNe. From the photometric 
aspect, however, SLSNe’s ®®Ni yields are usually rather 
small (e.g. Inserra et al. 201S) and the luminosities pro¬ 
vided by ®®Ni are outshone by magnetars or ejecta-CSM 
interaction, except for some PISN candidates that might 
be mainly powered by ^®Ni. Besides, ®®Ni yields are diffi¬ 
cult to be determined precisely owing to the lack of late¬ 
time spectra for most SNe. Hence, excluding a few PISN 
candidates, previous modelings for most LCs of SLSNe-I 
omitted the contributions of ®®Ni. 

There are some more complicated cases: for some 
SLSNe-H, e.g. SN 2006gy, adding several Mq of ®®Ni can 
make the fitting more reliable, the “^®Ni-|-interaction” 
model has been employed (Agnoletto et al. 2009; 
Chatzopoulos et al. 2012, 201S). No SLSN has involved 
with all three energy sources so far; Wang et al. (2015b) 
suggested that some luminous SNe Ic with peak magni¬ 
tudes — 20 mag must be explained by a unified sce¬ 
nario containing the contributions from ®®Ni and a mag¬ 
netar (“^®Ni-|-interaction” model) since they cannot be 


2 


powered solely by (see also Greiner et al. 2015) but 
the contribution from ®®Ni cannot be neglected (see also 
Bersten et al. 1967). 

Recently, an SLSN-I at redshift z = 0.3434, 

iPTF13ehe, might change the conclusions above. 
Yan et al. (2015) analyzed the light curve and spectra 
of iPTF13ehe and found that the peak bolometric lu¬ 
minosity of iPTF13ehe is ^ 1.3 x 10"^^ erg s“^. If this 
SLSN was powered by ®®Ni, more than (13 — 16)M0 of 
^®Ni must be synthesized to power the peak bolometric 
emission (Yan et al. 2015). However, the late-time spec¬ 
tral analysis indicated that the mass of ®®Ni synthesized 
is ~ 2.5 Mq. This amount of ®®Ni, although it cannot 
be ignored, is inadequate to power the peak bolometric 
luminosity. Therefore, Yan et al. (2015) argued that the 
early LC must be powered by multiple power sources. 

Furthermore, the presence of Ha emission lines in the 
late-time spectrum and the rebrightening in the late-time 
LC show an unambiguous signature for interactions be¬ 
tween the SN ejecta and a hydrogen-rich GSM. Hence, 
any model with single power source or double power 
sources is challenged in explaining the LC of iPTF13ehe. 
A triplc-power-source model should be seriously consid¬ 
ered. 

In this paper, we propose a triple-energy-source model 
for iPTF13ehe and try to understand its explosion mech¬ 
anism and LC. This paper is organized as follows. In Sec¬ 
tions 2 and 3, we fit the early-time and late-time LC using 
double and triple energy-source models, respectively. In 
Sections 4 and 5, we discuss and conclude our results. 

2. MODELING THE EARLY-TIME LC 

We start by studying the early LC and determining 
the power sources of iPTF13ehe. As pointed out by 
Yan et al. (2015), 2.5Mq of ®®Ni were synthesized 

while > 13 — I 6 M 0 of ^®Ni are required to account for the 
LC peak. Therefore, there should be other energy sources 
to aid the high peak-luminosity of iPTF13ehe. The early- 
time spectra shows no H emission lines, indicating that 
the ejecta-CSM interaction can be neglected in the early- 
time luminosity evolution. Hence, the energy source aid¬ 
ing the SN might be a magnetar spin-down input. The 
possibility that iPTF13ehe is a CCSN and can be aided 
by a magnetar was also pointed out by Yan et al. (2015), 
although they didn’t provide further analysis. Consider¬ 
ing the contribution from ®®Ni that cannot be neglected, 
the most reasonable model should contain contributions 
from both ®®Ni and a magnetar. Wang et al. (2015b) 
have constructed such a unified model for SNe Ic and 
explained three luminous SNe Ic. We here employ this 
model to fit the LC of iPTF13ehe. 

Supposing that the ratio of the r-band flux to bolomet¬ 
ric flux remains nearly constant around the peak and in 
the late-time, we scale the r-band luminosities to obtain 
the bolometric luminosities. The scaled bolometric lu¬ 
minosities are rough, but they can be regarded as an 
approximation of real values. These data are plotted in 
Fig. 1. In this section, we use some models to fit the 
early-time LC. The photospheric velocity and the opti¬ 
cal opacity k are fixed to ~ 13, 000 km s“^ (Yan et al. 
2015) and 0.2 cm“^ g“^, respectively. 

2 . 1 . Model 


The hrst of the models employed here is the ®®Ni- 
powered model. We find that when the ejecta mass 
Mej = 35 Mq, the gamma-ray opacity Kj ~ 0.018 cm“^ 
g“^, and the ®®Ni mass Mni ^ 19 Mq, the theoretical 
LC can matches the observational data at early times 
{t < 137 days), as can be seen in Fig. 1 (the dash-dotted 
line). The ^®Ni mass obtained here is slightly larger than 
the value ^ (13 — 16)M0 inferred by Yan et al. (2015). 

However, the ratio of Mm (~ 19 Mq) to Mej (35 Mq) 
is ~ 50%, significantly exceeding the upper limit of the 
values for CCSNe (the ratio of ^®Ni mass to the ejecta 
mass is < 20% (Umeda & Nomoto 2008)). Therefore, 
the ®®Ni-powered model can be excluded in explaining 
the data. 



Fig. 1.— Modeling iPTF13ehe’s early-time LC using ®®Ni, mag¬ 
netar, and hybrid (®®Ni -|- magnetar) models. Data are obtained 
from Yan et al. (2015). The early-time and late-time bolometric 
luminosities scaled from r-band ones are shown as the orange filled 
circles and orange unfilled circles, respectively. The dash-dotted 
line is the LC powered by 19 Mq of ®®Ni; the dashed line is the 
LC powered by a magnetar; the solid line is the LC powered by 
®®Ni and a magnetar. The short dashed line is the LC powered 
by 2.5 Mq of ®°Ni with full trapping of gamma-rays. The other 
model parameter values are shown in the text. 


2 . 2 . Magnetar Model and Magnetar+^^Ni Model 

In the magnetar-powered model, if the the initial rota¬ 
tional period of the magnetar Pq = 2.45 ms, the surface 
magnetic field B = 8 x 10^^ G, Kj ~ 0.018 cm“^ g“^, 
then the dashed line in Fig. 1 is plotted to obtain an 
equivalently good fitting. Since the spectral analysis has 
inferred that the SN explosion produced 2.5 Mq of 
®®Ni, we take into account the contribution from ™Ni 
and use the hybrid model, in which the initial period 
of the magnetar Pq is larger slightly, ~ 2.55 ms. From 
the fit, we can get other two important physical quan¬ 
tities relevant to the next section, i.e., the rise time 
A ^ 81 days and the kinetic energy of the SN explosion, 
Esn = ( 3 / 10 )Mejf ;2 ^ 3.5 x 10 ^^ 

2.3. Further Analysis for the Parameters of iPTFlSehe 

A caveat on the degeneracy of model parameters 
should be clarified here. The effective light-curve 
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timescale Tm = (2KMej/^t!c)^/^ (Arnett 1982) determines 
the width and therefore the rise time of the LCs of SNe, 
where c is the speed of light and ^ = 47r^/9 ~ 13.8 is 
a constant. If the value of nMe^/v is invariant, LCs 
around the peak are the same. Although the photo- 
spheric velocity v of iPTF13ehe is fixed to 13,000 km 
s“^, the degeneracy of k and A/ej is still difficult to be 
broken, because larger k requires smaller Mej for the 
same level of fittings. The optical opacity of H- and He- 
deficient ejecta has been assumed to be 0.06 cm^ g“^ 
(e.g., Valenti et al. 2011; Lyman et al. 2016), 0.07 
cm^ g“^ (e.g., Taddia et al. 2015; Wang et al. 2015b), 
0.08 cm^ g“^ (e.g., Arnett 1982; Mazzali et al. 2000), 
0.10 cm^ g“^ (e.g., Nugent et al. 2011; Inserra et al. 
2013; Wheeler et al. 2015), and 0.2 cm^ g“^ (e.g., 
Kasen & Bildsten 2010; Nicholl et al. 2014). We assume 
that K =0.2 cm^ g“^ throughout this paper and get 
the ejecta mass Mej ~ 35 Mq. If we assume a smaller 
K, e.g., K =0.1 cm^ g“^, then the ejecta mass becomes 
larger, Mej ~ 7OM0. The latter is consistent with the 
lower limit derived by Yan et al. (2015). Provided that 
iPTF13ehe is a genuine core-collapse SN, the ejecta with 
mass of 35Mq is preferred. 

3. MODELING THE LATE-TIME LC 

We have demonstrated that the hybrid model contain¬ 
ing the contributions from a magnetar and ®®Ni can well 
fit the early-time LC of iPTF13ehe, but the late-time re¬ 
brightening cannot be explained by such a hybrid model. 
Yan et al. (2015) found that the late-time spectrum of 
the SN displayed Ha emission lines which are usually 
present in SNe Iln and inferred that a Hydrogen shell 
ejected about 40 years before the SN explosion, locating 
at a distance of ~ 4 x 10^® cm from the SN progenitor. 
When the SN ejecta run into the shell, an interaction 
was triggered and thus the r-band LC was rebrightened 
and the Ha emission was prompted. Therefore, both 
the photometric and spectral features indicate that the 
interaction between the ejecta and CSM must play an 
important role in powering the late-time LC (Yan et al. 
2015). 

However, the exact time of the emergence of the ejecta- 
CSM interaction is uncertain due to the lack of the pho¬ 
tometric data at 60 — 170 days and the spectral data at 
13 — 251 days after the peak. Nevertheless, the fact that 
there is no rebrightening in ~ 56 days after the peak pro¬ 
vides a lower limit on the interaction onset time, while 
the rebrightening at ^ 180 days after the peak offers an 
upper limit. Therefore, we suppose that the rebrighten¬ 
ing started at ^ 56 — 180 days after the peak. Adopting 
A 81 days, we conclude that the interaction should be 
triggered between ^ 137 — 261 days after the SN explo¬ 
sion. Hence, it is reasonable to suppose that the onset 
time of interaction between the ejecta and the CSM-shell 
ti is 140 days. 

According to Equations (14)-(16) and (20) of 
Chatzopoulos et al. (2012), we reproduce the late-time 
LC powered by the ejecta-CSM interaction, seen in Fig. 
2. The fit parameters are the optical opacity of the H- 
shell k'= 0.34 cm^ g“^, the ejecta mass Mej = 35Mq, 
the SN kinetic energy Esn = 3.5 x 10^^ erg (these two 
parameters are derived from the previous section), the 
slope of the inner density profile of the ejecta <5 = 0, the 



Fig. 2. — Modeling iPTFlSehe’s whole LC by using the triple- 
energy-source (®®Ni + magnetar + interaction) model. Data 
are obtained from Yan et al. (2015). The bolometric luminosi¬ 
ties scaled from r-band ones are shown as the orange filled cir¬ 
cles. The dotted line is the LC powered by ®®Ni and a magne¬ 
tar; the dashed line is the LC powered by the ejecta-CSM interac¬ 
tion; the solid line is the LC powered by the triple-energy-source 
(“®®Ni-|-magnetar-|-ejecta-CSM interaction”) model. The values of 
other parameters are shown in the text. 


slope of the outer density profile of the ejecta n = 7, the 
power-law exponent for the CSM density profile s = 0 
(corresponding to a uniform-density shell), the density 
of the CSM pcSM = 10“^'^ g cm“^, the mass of the 
CSM Mcsm = 27M0, the onset time of interaction 
between the ejecta and the CSM-shell ti ^ 140 days, 
the conversion efficiency from kinetic energy to radiation 
e = 0.05 that is comparable to the typical value e ~ 0.1 
(Moriya et al. 2013b). 

The density of the CSM shell in our model is one order 
of magnitude larger than the value derived by Yan et al. 
(2015). Since the shock-ionized CSM shell is proba¬ 
bly optically thin to the Thomson scattering (Yan et al. 
2015), 'TThomson = k ' pw < 1, where w is the width of 
the shell, k'(Mcsm/ dvr i?^n;)r(; = k'(Mcsm/ dvr i?^) < 1, 
so Mcsm < AttR^/k'. Yan et al. (2015) supposed that 
the interaction time is 251 days after the peak, and de¬ 
rived that i? ~ 4 X 10^® cm. If w = O.IR, and k' = 0.34 
cm^ g“^, then Mcsm < 30Mq, n = p/rriH < 1/k'wtoh 
4 X 10® cm“® (Yan et al. 2015), and p < 0.74 x 10“^® g 
cm“®. 

This difference is partly caused by the uncertainty of 
the onset time of the ejecta-CSM interaction. We sup¬ 
pose that the interaction time is ^60 days after the peak 
and then have i? = ~ 1.6 x 10^® cm; if w = O.li?, 

and k'= 0.34 cm^ g“^, then Mcsm < 5Al0, n < 1.1 x 10® 
cm“®, and p < 1.84 x 10“^® g cm“®. As pointed out 
by Yan et al. (2015), the possibility that the CSM shell 
might be partly ionized rather than fully ionized can¬ 
not be excluded. If the ratio of ionized CSM to to¬ 
tal CSM is 1/5.4, then pcSM ~ 10“^"^ g cm“®, and 
AfcSM ~ 27M0, being consistent with the values of our 
fitting parameters. 

As shown in Fig. 2, the LC powered by the triple 
power-source model (“®®Ni-|-magnetar-|-ejecta-CSM in¬ 
teraction” model) is well consistent with the data. If 
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we ignore the contribution from ®®Ni, the LC can be re¬ 
produced by the magnetar and ejecta-CSM interaction. 

We caution that the bolometric LC is obtained by scal¬ 
ing the r-band LC and therefore is not accurate enough. 
The uncertainties in the bolometric luninosities prevent 
us from getting precise estimate of the model parame¬ 
ters. Nevertheless, this does not change the facts that the 
early-time LC is mainly powered by a magnetar rather 
than ^®Ni and that the late-time LC requires adding the 
ejacta-CSM interaction. More accurate bolometric cor¬ 
rections can pose more stringent constraints on the pa¬ 
rameters. 

4. DISCUSSIONS 

The ®®Ni production is usually ineffective after the 
explosion of a neutrino-driven CCSN. Radioactive ^®Ni 
in a neutrino-driven CCSN is mainly synthesized in 
two ways: a neutrino-driven wind and an explosive nu¬ 
clear burning of silicon (and other intermediate-mass el¬ 
ements) heated to a temperature exceeding 5 x 10® K. 
The neutrino-driven wind can produce \0~^Mq of 
®®Ni and can be neglected; the amount of ®®Ni synthe¬ 
sized in shock-heating silicon shell is related with the 
shock energy and therefore the kinetic energy of a CCSN. 
When the kinetic energy is ~ 10®^ erg, the mass of 
®®Ni is ~ (0.05 — O.1)M0. If the neutron star left be¬ 
hind the explosion is a millisecond magnetar, it can help 
the ejecta synthesize additional ®®Ni up to ^ 0.2Mq 
(Suwa & Tominaga 2015). The total amount of ®®Ni 
from these processes is < O.SMq, significantly less than 
2.5 Mq. 

Therefore, it seems that the model containing the con¬ 
tributions from 2.5 Mq of ®®Ni and a magnetar is prob¬ 
lematic. However, the possibility that the ®®Ni mass 
of iPTFlSehe has been overestimated cannot be com¬ 
pletely excluded. We can get some indirect clues from SN 
2007bi: Gal-Yam et al. (2009) proposed that the amount 
of ®®Ni synthesized by SN 2007bi is ^ 5M0 by matching 
the low-quality nebular-phase spectrum and concluded 
that SN 2007bi is a PISN, while Dessart et al. (2012) 
showed that the PISN model must produce cool photo¬ 
spheres, red spectra and narrow-line profiles which con¬ 
flict with observational properties of SN 2007bi, and con¬ 
cluded that SN 2007bi might be powered by a magnetar. 
This example indicates that the simple spectral analysis 
is not enough to determine the precise value of ®®Ni yield 
and a more detailed analysis is needed. In this paper, we 
can regard ~ 2.5Mq as the upper limit of ®®Ni yield of 
iPTFlSehe. 

We note that SN ejecta-CSM interaction might also 
be a promising scenario for explaining the early-time LC 
of iPTFlSehe. If the progenitor experiences some erup¬ 
tions and expels some discrete shells before its explosion, 
the ejecta can collide these shells at different epoches 
and the interactions between the ejecta and these shells 
can power the early-time LC. The intermittent spectral 
observations at early times might miss the spectral fea¬ 
tures. In this scenario, the putative millisecond magne¬ 
tar can be removed and the central remnant is therefore 
a black hole. However, how does the explosion synthe¬ 
size ~ 2.5Mq of ®®Ni in this scenario remains an open 
problem. 

5. CONCLUSIONS 


We have analyzed and htted the LC for iPTFlSehe 
which is an SLSN I with peak bolometric luminosity 
~ 1.3 X 10^^ erg s“^. Using the ^®Ni-powered model, 
we found that ^ 19M0 of ®®Ni are required, confirming 
the conclusion reached by Yan et al. (2015). Spectral 
analysis performed by Yan et al. (2015) shows that the 
SLSN synthesized ~ 2.5Mq of ®®Ni but this amount of 
®®Ni is inadequate to power the peak luminosity. 

We proposed that the most reasonable model account¬ 
ing for the early-time LC is the magnetar-dominated 
model containing the contributions from a magnetar with 
Po = 2.45 - 2.55 ms and H = 8 X 10^^ G and < 2.5Mq 
of ®®Ni. The LCs reproduced by this model are in good 
agreement with early-time observations of iPTFlSehe. 
The energy released by the magnetar is dominant in pow¬ 
ering the optical LC of iPTFlSehe. 

While proposing this hybrid model containing triple 
energy-sources, we also note that this model is not unique 
model of explaining iPTFlSehe since the early LC might 
also be powered by the ejecta-CSM interaction and ®®Ni. 
Since the CCSNe usually produce < O.3M0 of ®®Ni, ~ 
2.5Mq of ®®Ni is problematic in both “magnetar+®®Ni” 
model and “interaction-|-®®Ni”. We regard 2.5Mq as 
the upper limit and call for more detailed analysis. 

An important spectral feature owned by iPTFlSehe 
is Ha emission lines in its late-time spectrum. These 
lines, together with the rebrightening of the late-time 
LC (Yan et al. 2015), indicate that the interaction be¬ 
tween the SN ejecta and the CSM has been triggered. 
The hybrid model involving ®®Ni and magnetar cannot 
account for the late-time rebrightening. We hnd that the 
ejecta-CSM interaction model can explain the rebright¬ 
ening and is the natural requirement accounting for the 
Ha emission lines. Therefore, we use the model contain¬ 
ing a magnetar, ®®Ni, and ejecta-CSM interaction to fit 
the whole LC. 

Finally, we discuss the explosion nature of iPTFlSehe. 
The physical properties of the progenitors, the explosion 
mechanisms and energy-source mechanisms of SLSNe 
are still ambiguous and controversial (see Dessart et al. 
2012; Nicholl et al. 2013, for further analysis). As em¬ 
phasized by Yan et al. (2015), the massive Hydrogen- 
rich shell must have been expelled by the “pulsational 
pair-instability (PPI)” mechanism (Heger et al. 2003; 
Woosley et al. 2007; Pastorello et al. 2008; Chugai 2009; 
Chatzopoulos & Wheeler 2012). The ejecta colliding 
with CSM shell contains < 2.5Mq of ®®Ni and > 3.5 x 
1052 gj.g qJ kinetic energy. This ejecta is not a shell 
following the hydrogen-rich shell, because Heger et al. 
(2003) pointed out that the shell expelled by PPI con¬ 
tains no ®®Ni and its kinetic energy is up to several times 
10®^ erg. The kinetic energy of the ejecta of iPTFlSehe is 
at least one order of magnitude higher than that of each 
PPI pulsed shell. Therefore, we propose that iPTFlSehe 
itself is a genuine CCSN rather than a PPISN candidate, 
while the hydrogen-rich massive shell could be caused 
by a PPI pulse ~ 20 years before the CCSN explosion 
{St ~ vti/vcsM ^ 13,000 km/s x 140 days/300 km/s ^ 
17 years, the velocity of the CSM shell unsM was inferred 
by Yan et al. (2015)). 

Yan et al. (2015) estimated that at least 15% of 
SLSNe-I may own such interaction features for late-time 
spectra. We expect that future studies focusing on sim- 
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ilar SLSNe should provide further insight into their ex¬ 
plosion and energy-source mechanisms. 
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